Role of boundary constraints in DNA cyclization 
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We modify the classical Shimada-Yamakawa theory of DNA looping by generalizing the form of 
boundary constraints. This generalization is important in the context of DNA cyclization exper- 
iments since it mimics the reduced local rigidity of the "nicked" DNA loop. Our results indicate 
that the non-trivial boundary constraints may be responsible for the existing dramatic discrepancy 
between various DNA cyclization experiments. The developed effective Hamiltonian method may 
be extended to even broader class of DNA looping problems. 
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Loop formation in double-stranded DNA (dsDNA) is 
essential for such important biological processes as reg- 
ulation of gene expression and DNA packaging into nu- 
cleosomes. In the first case, the loop is induced by the 
interactions between transcription factor proteins bound 
to different sites along the DNA chain [H-@]. In the case 
of nucleosome, DNA wraps around the near-cylindrical 
histone octamer. The classical theory of looping, based 
on the elastic description of DNA, was proposed more 
than two decades ago by Shimada and Yamakawa (SY) 
[3| ■ There are however multiple indications that the orig- 
inal theory is not sufficiently adequate for describing the 
real experimental situation. 

Partially, the deviations can be attributed to the com- 
plexity of the actual in vivo problem. However, a large 
discrepancy is also reported in recent in vitro experiments 
on DNA cyclization [f| @ . In these experiment, looping 
is induced by hybridization of mutually complementary 
ssDNA ends of the chain (see Figure [T]). Even in this rel- 
atively simple case, the agreement between theory and 
experiment is a highly controversial issue. 

On the one hand, older cyclization experiments by 
Shore et al do agree with SY model 7][8]. In fact, those 
data were used in the original SY paper to support their 
model. On the other hand, Cloutier and Widom reported 
that the discrepancy between experimental looping prob- 
ability and SY results may reach two to three orders of 
magnitude [5|. Their paper was followed by the work 
of Vologodskii lab, in which the agreement between the- 
ory and experiment was confirmed [9(. More recently, 
Cloutier and Widom reiterated their claim, and also re- 
ported that the twist-related oscillations of the looping 
probability are strongly reduced compared to the SY pre- 
diction 0. This controversy inspired a new interest to 
the problem among theorists. In particular, it was sug- 
gested that strong bending may induce local structural 
defects such as "kinks" or ssDNA "bubbles" (which act 
as "soft" kinks) [10|-[l^|. At present, there is no direct ev- 
idence to support any of these models, and they do not 
resolve the conflict between the different experiments. 

In this communication, we propose an alternative ex- 
planation to the existing controversy. We argue that the 
DNA cyclization involve more complicated boundary (i.e. 
terminal) constraints then it is traditionally believed. In 




FIG. 1: Scheme of DNA cyclization experiment. Looping is 
caused by hibridization of mutually complementary ssDNA 
ends. 



order to solve the problem with the modified boundary 
conditions, we develop an effective Hamiltonian method. 
The approach is potentially applicable beyond the scope 
of this work, e.g. for study of protein-mediated looping 

The SY theory of looping starts with dsDNA mod- 
elled as an elastic rod subjecxted to thermal fluctua- 
tions, which is essensially equivalient to Worm Like Chain 
(WLC) model 
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Here, the first term describes the bending elasticity 
with modulus proportional to dsDNA persistence length 
l p w 50nm, the second term represents torsional elastic- 
ity, with the modulus proportional to twist persistence 
length, l t ~ lOOnm. t is unit tangent vector of the chain, 
and ip is twist angle, both functions of position along the 
chain, s. 

From the experimental and biological points of view, 
the most important characteristics of the looping prob- 
lem is the so called J- factor. It has a meaning of an effec- 
tive concentration of one end of the loop in the vicinity of 
the other, at the open configuration. J-factor can be re- 
lated to equilibrium probability for the loop to be closed, 
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Here e is the mutual affinity of the loop terminals, and 
Co = 1M is the standard reference concentration. Affin- 
ity e can be independently determined from the experi- 
ment on dimerization of free reacting terminal groups in 
solution. 

Shimada and Yamakawa calculated J-factors for two 
important cases: circular loops with completely aligned 
ends (t(L) = t(0)): 
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and loops with unconstrained orientations of the end seg- 
ments: 
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These results are obtained in the limit of short loops. 
However, their range of applicability extends up to L ~ 
10lp (i.e. 1500 bp). The effects of torsional constraints 
which are not included in Eqs.([3])-([4]), result in an addi- 
tional factor with a pronounced oscillatory behavior with 
the period of DNA helical turn, 10.5 bp. 

DNA loop formed in a reversible cyclization experi- 
ments (i.e. before ligation) is not identical to a circular 
DNA. It is "nicked" in two points corresponding to the 
ends of the DNA strands (see Figure [T]). These singular 
points are expected to have a greater flexibility then the 
rest of the chain [l(|. The effective local flexibility must 
strongly depend on the base-stacking interactions. We 
will describe the coupling between orientations of the end 
segments of the loop with the following minimal model: 
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Here 9 is the angle between directions of tangent vectors, 
t (0) and t (L), and Aip is the relative twist of the two 
segments. The two parameters, k and k' should be of 
the same order for the given sequence. Since the char- 
acteristic stacking energy is of order of kT, we expect 
k ~ k' ~ 1. 

In order to calculate J— factor with this modified 
boundary constraints, we develop an effective Hamilto- 
nian method (EHM) for the problem. Effective Hamil- 
tonian can be introduced as a free energy of the loop 
parameterized with orientations of its end segments. If 
we neglect any torsional constraints (i.e. assume k' = 0), 
the Hamiltonian of a closed DNA loop can be written 
as a function of azimuthal and polar angles, 9 and ip of 
tangent vector t (L), with respect to the vector triad at 
the other end of the loop, s = 0: 
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In order to construct Hi oop (6,<p), we first discuss the 
ground state energy of the loop. Since the original 
Hamiltonian of WLC model is formally equivalent to La- 
grangian of symmetrical spinning top, finding its ground 
state is an integrable mechanical model. However, find- 
ing the exact solution for a particular set of boundary 
conditions, requires an inversion of incomplete Elliptical 
functions, and therefore it is not practical [13| . Instead, 
we consider the problem in the vicinity of circular loop 
configuration (which corresponds to t (0) = t(L)). In 
this limit, the ground state energy of a 2D (planar) loop 
can be written as an expansion in 9: 
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where the exact values of the coefficients are: (3 = 2n; 
7 = 3. Already this expression gives an excellent de- 
scription of the global behavior of function Ei oop (9), not 
limited to the near vicinity of circular loop (9 = 0). In 
particular, it predicts the minimal energy, E m - ln /kT = 
(47T 2 /3)Z p / L ps 13.16^/L, which is very close to the 
results of complete numerical solution of the problem: 
Emin/kT w 14.04l p /L. Nevertheless, due to the strong 
exponential dependence of J-factor on the elastic energy, 
we need a further improvement of the analytic formula 
for Eioop (9)- We achieve this by treating the coefficient 7 
in as a free parameter, and adjusting it to value 7 = 3.46. 
This leads to an exact matching of the minimum energy. 
The overall behavior of approximate expression, Eq.([7]), 
becomes nearly indistinguishable from the exact result 
(see Figure [5]). 




FIG. 2: Comparison of the Effective Hamiltonian (solid line), 
to the exact result for the elastic energy of a loop (squares). 
Dashed line shows the extrapolation from point 8 = 0, with- 
out any adgustment to model parameter 7 (i.e 7 = 3). 

Since the change in the elastic energy is the dominant 
correction to the loop free energy, we can write the effec- 
tive Hamiltonian of the loop in the following form: 
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Here, the earlier expression, Eq. Q, has been general- 
ized for the case of 3D loops, and additional fine-tuning 
parameter, 7' has been introduced. Fq is free energy of 
the circular loop. 

The combined effective Hamiltonian Eq. ©, reaches 
minimum at 8 = 8* = — 2t:/ (7 + nL/l p ), and tp = 0. 
We now expand it in the vicinity of that point: 
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This quadratic expansion allows one to calculate the over- 
all free energy of the looped state in Gaussian approxi- 
mation, and hence obtain J-factor for arbitrary coupling 
n : 
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The result of this calculation can be well approximated 
by the following analytic expression (the torsional effects 
are omitted): 
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This form generalizes of the original SY result The pa- 
rameter 7' was tuned to the value 7' = —7/2, to achieve 
a nearly exact matching with both SY limits, Eqs. © 
and (4}. 

To complete our calculation we now need to include 
effects of the torsional constrain. For short loops, the 
corresponding effective Hamiltonian is the sum of the 
elastic energy of the twisted chain and the torsional part 
of H en d'- 
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Here ./V is an integer linking number, and h is helix repeat 
of dsDNA. This leads to an additional torsion-related 
factor in the final result: 
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FIG. 3: Comparison of the theoretical result, Eq. Q with ex- 
perimental data from Refs. Q| (diamonds) and (squares). 
The two curves correspond to weak (« = 1, k' — 2) and strong 
(k = k' = 20) orientational couplings, respectively. 



The calculated J-factor is shown in Figure[3l as a func- 
tion of the loop size. Remarkably, the two conflicting 
sets of experiments are both consistent with the model. 
In particular Cloutier-Widom and Vologoidskii lab data 
correspond to the regimes of weak (k = 1, n' — 2) and 
strong (k = k' = 20) terminal coupling, respectively. 
As expected, the two coupling parameters, k and k! are 
strongly correlated. This variation in value of k is reason- 
able since the effective rigidity of a DNA "nick" is likely 
to have an exponential dependence on the local stacking 
energy. This strong dependence of the cyclization proba- 
bility on the local sequence should not be confused with 
another sequence dependent effect associated with inho- 
mogeneous intrinsic curvature and bending modulus of 
the chain In order to separate the two effects, the 
sequence of the mutually complementary terminal groups 
and that of the rest of the chain should be varied inde- 
pendently in the future experiments. 



One can make several important conclusions based on 
our results. First, the traditional "circular loop" mod- 
elling of DNA cyclization is only justified for very strong 
coupling, k > 10. Second, in the regime of moderate 
coupling (k ~ k' ~ 1) the effects of the two constraints 
are rather different. The overall shape of J-factor curve 
follows that of an orientationally unconstrained loop, Eq. 
(2]), but the torsion-related oscillations, although greatly 
reduced, remain rather prominent. This is indeed consis- 
tent with the results of Ref . [6| . 
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